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The low-temperature shrinking of the vortex core (Kramer-Pesch effect) is 
studied for an isolated single vortex for chiral p-wave and s-wave supercon- 
ducting phases. The effect of nonmagnetic impurities on the vortex core ra- 
dius is numerically investigated in the Born limit by means of a quasiclassical 
approach. It is shown that in the chiral p-wave phase the Kramer-Pesch effect 
displays a certain robustness against impurities owing to a specific quantum 
effect, while the s-wave phase reacts more sensitively to impurity scattering. 
This suggests chiral p-wave superconductors as promising candidates for the 
experimental observation of the Kramer-Pesch effect. 

PACS numbers: 74.25. Op, 74.20. Rp, 74.70. Pq 



1. INTRODUCTION 

Much attention has been focused on vortices in nature^ especially for 
the quantized vortex in fermionic superfluid and superconducting systems) 2 * 3 * 4 * 
One of the fundamental physical quantities of the quantized vortex is the 
radius of vortex core. Kramer and PesclP have pointed out theoretically 
that the radius of vortex core decreases proportionally to the temperature T 
at low temperatures, much stronger than anticipated from the temperature 
dependence of the coherence length. This anomalously strong shrinking of 
the vortex core, the so-called Kramer-Pesc h (KP) effect,* 5 - *^* occurs when the 
fermionic spectrum of vortex bound statej 7 -*^* 9 * 10 * 11 * 12 ^ crosses the Fermi 
levelPl The temperature dependence of the vortex core has been theoreti- 
cally investigated in the case of superconductors! 9 * 15 * 16 * 17 * 18 * 19 * 20 * 21 * 22 * 23 * The 
low-temperature limit of the vortex core radius was discussed also for dilute 
Fermi superfluids^ and superfluid neutron star matter.^ 

There are several length scales which characterize the radius of vortex 
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core! I I I One of them is the coherence length £(T). The pair potential 
A(r) depressed inside the vortex core is restored at a distance r ~ £(T) 
away from the vortex center r = 0. 27 However, £(T) is almost temperature 
independent at low temperatures. Another length scale is related to the 
slope of the pair potential at the vortex center, which is defined as 

T = t~7 r lim — — . (1) 

£1 A(r —> oo) r^o r 

The KP effect means £i(T) oc T for T — > 0, while the pair potential A(r) is 
restored at a distance r ~ £ (3> £1) at low temperatures.^ Since the spatial 
profiles of the pair potential A(r) and the supercurrent density j(r) in the 
vicinity of the vortex center are related with each other through the low- 
energy vortex bound states, the length g i scal es with the distance r = ro at 
which \j(r)\ reaches its maximum value ! 5 ' 2b ' 28 ' Therefore, the KP effect gives 
rise to ro ~ £1 — > (T — ► 0) linear in T. 

Recently, /iSR experiments we re perform ed on the s-wave superconduc- 
tor NbSe2 to observe the KP effect P^* 2 - 8 ' 29 ' 30 ' The experimental data of spin 
precessions, which correspond to the Fourier transformation of the Redfield 
pattern (the magnetic field distribution) of a vortex lattice, are fitted by 
a theoretical formula,^ extracting the information on the spatial profile of 
the supercurrent density j(r) around vortices) 26 * 28 * While a shrinking vortex 
core was observed, it was weaker than the theoretical expectation of the KP 
effect for a clean superconductor. That is, indeed the observed vortex core 
radius ro seemingly shrank linearly in T, but was extrapolated towards a 
finite value in zero-temperature limit, indicating a saturation of the KP core 
shrinkage at a low temperature. The measurement of the vortex core radius 
ro by /iSR was also performed for CeRu2 with a similar result P^- Since the 
KP effect is directly connected with the low-energy vortex bound states, the 
energy level broadening due to impurities may give rise to a modification 
of the KP effect,^ in particular, a saturation as observed. Impurities ex- 
ist inevitably in any solid state material, so that such a saturation effect is 
not unlikely to occur. Even it turns out that a rather small concentration 
of impurities, leaving the material moderately clean, would have a strong 
influence as observed in experiments!^* 2 - 8 * 29 * 30 * 32 ' 

There are several factors which influence the behavior of vortex core 
radius, (i) Impurity effects^ (ii) The discreteness of the energy levels of the 
vortex bound states! 5 ' 2 " ' 21 ' 22 ' (iii) Vortex lattice effectsP^EHESEiElESEzI ( iv ) 
Fermi liquid effects Ff. 3S (v) Antiferromagnetic correlations induced inside 
the vortex core as suggested for cuprate superconductors!*^ (vi) The presence 
of multiple gaps in multi-band superconductors! 4 ^ ' 41 ' 

In this paper, we investigate the temperature dependence of the vortex 
core radius £i(T) incorporating the effect of nonmagnetic impurities on the 
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level of the Born approximation, for both a chiral p-wave and an s-wave 
superconductor. For this purpose, we set up a model of single vortex in 
a two-dimensional superconductor, where the chiral p-wave pairing has the 
forrrP^EHI d(k) = z(k x ± ik y ) as, for example, realized in Sr2Ru04. 

The examination of the temperature dependence of £i(T) suggests that 
under certain conditions the chiral p-wave state shows a more robust KP ef- 
fect against impurities than an s-wave state. This behavior is connected with 
the compensation of the intrinsic phase structure and the vortex phase wind- 
ing, if chirality and phase winding are oppositely oriented! 2 ^ * 44 * 45 * 4 ^ 4 ^ 4 ^ 4 ^ 5 ^ * 51 * 52 * 
Therefore, we argue that chiral p-wave superconductors might be better can- 
didates for the experimental observation of the KP effect. The chiral p-wave 
superconductivit y has been proposed for Sr2Ru04 with rather stron g exper- 
imental evidence? 42 * 4 ^ and more recently for Na^CcC^ ■ t/H^OP^^ 

The paper is organized as follows. In Sec. 2, the self-consistent system 
of equations for the superconducting order parameter and the impurity self 
energy is formulated on the basis of the quasiclassical theory of supercon- 
ductivity. In Sec. 3, the systems of the s-wave vortex and the chiral p-wave 
vortex are described. The numerical results for the vortex core radius 
are shown in Sec. 4. The summary is given in Sec. 5. 



2. QUASICLASSICAL THEORY 

Our theoretical analysis of the KP effect will be based on the quasiclas- 
sical theory of superconductivity,* 55 * 56 * which allows us to take inhomogeneous 
structures as a vortex into account and at the same time to deal with impu- 
rity scattering on a straightforward way. 

We start with the quasiclassical Green function in the presence of non- 
magnetic impurities, 

g{iu m r, k) = -in ( _ 9 . f] % J \ , (2) 

which is the solution of the Eilenberger equation 

zv F (k) • X7g + [iu) n T 3 - A - a, g] = 0. (3) 
Here, A is the superconducting order parameter 

and a denotes the self energy correction due to impurity scattering 

a(iuj n ,r,k) = ( 011 0x2 ) . (5) 
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The Eilenberger equation is supplemented by the normalization condition 
g(iu> n ,r, k) 2 = — 7r 2 iP^'The vector r is the real-space coordinates and the 
unit vector k represents the sense of the wave number on the Fermi surface. 
vp(k) is the Fermi velocity, u) n = irT(2n + 1) is the fermionic Matsubara 
frequency (with the temperature T and the integer n), T3 is the third Pauli 
matrix in the 2x2 particle-hole space, and the commutator [a, b] = ab — ba. 
We will consider an isolated single vortex in extreme type-II superconductors 
(Ginzburg-Landau parameter k 3> 1), and therefore the vector potential is 
neglected in the Eilenberger equation (j3J). Throughout the paper, vectors 
with the upper bar denote unit vectors and we use units in which h = k& = 1. 
We now define an alternative impurity self energy S as 

VV- f\ ( S d £12 \ ( lO^ll — <J 22) Cl2 \ , R \ 

£zw n ,r,k = = 2v I 6 ) 

The original impurity self energy @ can be expressed as 

a = ± + <fl + !fi. (7) 

Hence, we rewrite the Eilenberger equation (J3J) as 

ivp(k) • Vg + [iQ n T 3 -A,g} = 0, (8) 
with the renormalized Matsubara frequency 

iuj n = iu) n - S d , (9) 
and the renormalized superconducting order parameter 

MV-*o A+ » Sl2 )- 

We restrict here to s-wave scattering at the impurities. The single- 
impurity t-matri^^ is then calculated as 

i(iu n ,r) =vl + N v(g(iuj n , r, k)ji(iuj n , r), (11) 

where v is the impurity potential for the s-wave scattering channel, iVo is 
the normal-state density of states at the Fermi level, and the brackets (• • •) 
denote the average over the Fermi surface. The impurity self energy a is 
given by 

a Tt' V 

a(iuj n ,r) = mt(iu n ,r) = — 



1 + N v(g(iu n ,T,b)) , (12) 
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where the denominator is 

D=l + (irN v) 2 [(g) 2 + {f)(f% (13) 

and n\ is the density of impurities. The scattering phase shift 5q is defined 
by tanJo = —ttNqv. In this paper, we investigate the Born limit (<5o <C 1). 
The impurity self energy (|12|) in this limit becomes 



a(iu n ,r) = mvi + —(g(iu) n ,r,k) 

- "*I + r.(r$ W), (14) 

where we have defined the impurity scattering rate in the normal state as 
T n = l/2r n = irriiNov 2 . The mean free path I is defined by I = VpT n = 
VF/2T n . From Eqs. © and (fill) , we obtain the self-consistency equations 
for S as 



T, d (iui n ,r) = -iT^giiu^rik)), (15) 
Ei 2 (iu n ,r) = r n (/(io; n ,r,k)), (16) 



E 2 i(w n ,r) = -r n |/ f (ia; n ,r,k)|. (17) 
The self-consistency equation for A, called gap equation, is given as 

A(r,k)=vrT 5 F(k) £ (**(£)/(««;„, r, k')Y (18) 



— U) c <Ul n <U) c 



where the cutoff energy is u c , the pairing interaction is defined as gF(k)F* (k') 
with the coupling constant g given by 

We define T c q as the superconducting critical temperature in the absence 
of impurities. Finally, the system of equations to be solved consists of the 
Eilenberger equation (jHJ) and the self-consistency equations for the impurity 
self-energies [Eqs. (|15|) - (|17|) ] and for the superconducting order parameter 
[Eq. CHD]. 
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3. 5- WAVE AND CHIRAL P-WAVE VORTEX SYSTEMS 

In this study, the system is assumed to be an isotropic two-dimensional 
conduction layer perpendicular to the vorticity along the z axis. In the 
circular coordinate system we use r = (r cos <fi, r sin <fi) and k = (cos 8, sin 6). 
We assume a circular Fermi surface and vp(k) = i>Fk = (vf cos 6, vp sin 6) in 
the Eilenberger equations (j3J) and 0. The average over the Fermi surface 
reads: (•••) = ■ ■ • cZ0/2tt. 



3.1. Pair Potential 

The single vortex is situated at the origin r = 0. In the case of the s- 
wave pairing, the pair potential (i.e., the superconducting order parameter) 
around the vortex is expressed as 

A s (r) = A s (r)e^, (20) 

where we can take A s (r) in the right-hand side to be real because of axial 
symmetry of the system. The s-wave pairing means P(k) = 1 and the gap 
equation (JT8|) reads: 

A s (r)=7rTg s £ / —f( iu ; n ,r,9), (21) 

-u c <w„<ta c " 

where g s follows Eq. (|19j). 

On the other hand, in the case of the chiral p-wave pairing, -F(k) = k x ± 
iky = exp(±i#). The pair potential around the vortex has two possible forms 
depending on whether the chirality and vorticity are parallel or antiparallel 
each other pD 51 l 52 l Thus, there are two kinds of vortex. One form is 

A( n )(r,fl) = A^(r)e +W + A^\r)e- id 

= APWe^'+A^We'W, (22) 

where the chirality (related to the phase "+#") and vorticity ("—0") are 
antiparallel ("negative vortex"). The other form is 

A( p )(r,0) = A ( f\r)e +W + A ip \r)e- ie 

= AfWe^+A^We 4 ^, (23) 

where the chirality ("+#") and vorticity ("+(/ ) ") are parallel ("positive vor- 
tex"). We have assumed that A^ n ' P ' ) is the dominant component in both Eqs. 
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(|22[) and (|23|) . and the other component A_ is the minor one induced with 
smaller amplitudes inside the vortex core. The axial symmetry allows us to 
take A±' p \r) to be real in each second line of Eqs. ()22j) and (|2"5|) . Far away 
from the vortex, the dominant component A+' p \r — > oo) is finite and the 
induced minor one A^' p ^(r — > oo) — > 0, namely 

A ( n )( r oo,<&0) = Aj°(r -> oo)e^~^, (24) 



A(p)( r 00,0; 0) = A^ p) (r -> oo)e i(e+< «. (25) 

In the clean limit, Af* p) (r -» oo) = A B csCO [A BC s(r) is the BCS gap 
amplitude]. The gap equation (fl8|) reads now: 

Ai n ' p) (r)=vrT 5p £ J^^o f{l ^ ^ (26) 
where g p follows Eq. (fTT))) . 



3.2. Axial Symmetry and Boundary Condition 

The numerical calculation of the self-consistent pair potential A and 
impurity self energy E requires to restrict ourselves to a finite spatial region, 
which we choose axial symmetric with a cutoff radius r c . Therefore, it is 
necessary to fix the values of A and E for r > r c outside the boundary 
r = r c when solving the Eilenberger equation ©. 

We set, for the pair potential (|2Ujl of the s-wave vortex, 

As(r>r C) 0) = A B (r = r c )e^. (27) 
For the pair potential (|22|) of the chiral negative p-wave vortex, we set 

A^(r >r c ,<t>;Q) = A^\r = r c )e i( - e ~^ , (28) 
and for the pair potential (|23|) of the chiral positive p-wave vortex 

AW (r > r c , 6) = A^ p) (r = r c )e i(9+ ^ , (29) 

while AL n ' p) (r > r c ) = 0. 

Next we consider the symmetry property and boundary condition of the 
impurity self energy E. In the Eilenberger equation (jHJ), £ appears in Eqs. 
@ and (fTU|) in the form of iio n — Ed, A + E12, and A* — E21. Owing to the 
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axial symmetry of our model, an axial rotation of the system {(ft — > (ft + a 
and 9 — > 9 + a) leads to a transformation of the impurity self energies Ed, 
E12, and E21 in the same manner as iu n , A, and A*, respectively. The 
Matsubara frequency iuj n is invariant under the axial rotation, and therefore 
Ed(iu; n , r, (ft) has no azimuthal (ft dependence. Hence, we set 

E d (iLJ„, r, (ft) = E d (iw„, r), (30) 

T, d (iu n , r > r c ) = E d (icj n , r = r c ), (31) 

both in the case of the s-wave vortex and the chiral p-wave vortex. 

For the s-wave vortex, the pair potential A s (r) in Eq. (|2()jl transforms 
under an axial rotation as A s (r, cj)+a) = A s (r) exp[i(0+a)] = A s (r, (ft) exp(ia), 
and A*(r, + = A*(r, 0) exp(— ia). This rotation means that A s + Ei2 — > 
(A s + Ei2)exp(ia) and A* — E21 — > (A* — E21) exp(— ia). Thus, the off- 
diagonal impurity self energies Ex2,2i(i^n) r , 4>) have to possess an azimuthal 
<f> dependence as 

Ei 2 (iuj n ,r,<l)) = Ei 2 (iw„,r)exp(i0), (32) 

E 2 i(iw n ,r,0) = E 2 i(iwn, , r)exp(-z0). (33) 
We set, for the s-wave vortex, the boundary condition as 

Ei 2 (iw„, r > r c ) = Z 12 (iuJn, r = r c ), (34) 
E 2 i(iw„,r > r c ) = E 2 i(iu; ri ,r = r c ), (35) 

because far away from the vortex core the anomalous Green functions aver- 
aged over the Fermi surface, which appear in Eqs. (|16|) and 1|17|) . are gener- 
ally nonzero and their amplitudes are spatially uniform, owing to the s-wave 
pairing symmetry.^ 

For the chiral negative p-wave vortex, the pair potential A^ has the 
following symmetry properties. For an axial rotation, the pair potential 
A( n )(r, 9) in Eq. (J22J) transforms as A( n )(r, (ft + a; 9 + a) = A (n \r, (ft; 9) and 
A( n )*(r, (ft + a; 9 + a) = A^*(r, (ft; 9), namely invariant. Therefore, we find 
that AW + E12 -» AW + E12 and A^* - E 2 i -> A( n )* - E 2 i and the 
off-diagonal impurity self energies E 12, 21 (i^n, f, <ft) are n °t ^-dependent: 

Y,i 2 (iu) n , r, (ft) = Ei 2 (m> n , r), (36) 

Y, 2 i{iu) n ,r,(j)) = E 2 i(io; n ,r). (37) 
We set, for the chiral negative p-wave vortex, the boundary condition as 



Ei 2 (iw n ,r > r c ) = 0, 



(38) 
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T, 2 i{iu; n ,r > r c ) = 0, (39) 

because far away from the vortex core the anomalous Green functions aver- 
aged over the Fermi surface are zero, owing to thep-wave pairing symmetry.^ 
On the other hand, the chiral positive p-wave vortex behaves differently 
under axial rotation. The pair potential A^ p \r,6) in Eq. (|23[) transforms 
as + a; 9 + a) = A(P)(r,0;0)exp(2ia) and A<p> (r, 4> + a; + a) = 

A( p )*(r,c/>;/9)exp(-2ia), such that A^) + E u -» (AW + E12) exp(2ia) and 
A(p)* - S21 (A( p )* - £21) exp(— 2ia). Here, the off-diagonal impurity self 
energies £12,21 («^n> f-> 4>) depends on <p like 

Y, 12 (iuj n , r, 4>) = T, 12 (iu n , r) exp(2i<j>), (40) 

£ 2 i(iw n , 0) = E 2 i(jw n) r) exp(-2i0). (41) 
We set, for the chiral positive p-wave vortex, the boundary condition as 

Si 2 (icj„,r > r c ) = 0, (42) 

E 2 i(m>n,r > r c ) = 0, (43) 
as in the case of the negative vortex. 



4. VORTEX CORE RADIUS £i(T) 

The vortex core radius £i(T) defined in Eq. Q is obtained from the 
spatial profile of the pair potential calculated self-consistently. We solve the 
Eilenberger equation (JHJ) by a method of the Riccati parametrization j lu * 58 land 
then iterate the calculation, until self-consistency is reached, with the self- 
consistency equations of the impurity self energies [Eqs. I|15 |) -(|17j ) ] and that 
of the pair potential for the s-wave vortex [Eq. (|21|) ] or for the chiral p-wave 
vortices [Eq. (|26|)]. We use an acceleration method for iterative calculations 
to obtain sufficient accuracy.^ When solving the Eilenberger equation 0, 
we use the boundary conditions for the pair potential and the impurity self 
energies described in Sec. 3.2. 

We show now the results of £i(T) for the three cases introduced above. 
We set the cutoff energy uj c = 10 Ao and the cutoff length r c = 10£o- Here, 
Ao is the BCS gap amplitude at zero temperature in the absence of impu- 
rities, and £0 = ^f/Ao- The lowest value of the temperature at which £1 is 
computed is T = 0.02T c o. 

We have checked an influence of the finite system size, comparing results 
of £1 obtained for two cutoff lengths r c = 10£o and 20£o in the case of the 
clean limit. As a result, it was found that those results well coincide each 
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other below T ~ 0.7T c o (the deviations in £1 are less than 0.003£o and 
negligibly small), while slight deviations appear above T ~ 0.8T c o, but are 
less than 0.014£o (T < 0.9T c o) and almost invisible for plots with the same 
plot scale as in the following figures. 



4.1. .S-Wave Vortex 

In Fig. 1, we show the vortex core radius £i(T) for the s-wave vortex as 
a function of the temperature T for several values of the impurity scattering 
rate r n . The critical temperature T c remains unaffected by the nonmagnetic 
impurities, namely T c = T c q for any values of T n . 

At high temperatures, we see in Fig. 1 that the vortex core radius £i de- 
creases with the increase of the impurity scattering rate r n . This is because 
the coherence length £, over which the pair potential significantly changes, 
shrinks with the decrease of the quasiparticle mean free path (oc l/r n )P3 
The coherence length £ is a distance at which the pair potential is restored 
far away from the vortex center. The vortex core radius £i defined in the 
vicinity of the vortex center is dominated by this r n dependence of the co- 
herence length £ in this temperature regime. 

Pronounced T n dependence appears also at low temperatures. For r n — 
(the clean limit), the vortex core radius £i decreases linearly in T, as 
expected for the KP effect. In Fig. 1(a), we also show £i(T) for finite values 
of the impurity scattering rate, r n = 0.1 Ao and Ao- At low temperatures, 
the vortex core radius £i increases with the increase of r n , in contrast to 
the high-temperature behavior. This increase of £i indicates the saturation 
feature of the KP effect due to impurities P The low-temperature vortex core 
radius £i expands by introducing impurity scattering but still remains much 
smaller than £. For relatively small r n (< Ao), the decrease of the coherence 
length £ mentioned above has little influence on this expansion of the vortex 
core radius £i (<S £). For larger r n towards the dirty limit, however, the 
decrease of the coherence length £ begins to influence the vortex core radius 
£i, and £i begins to decrease with growing r n as seen in Fig. 1(b). 

In Fig. 1(b), the solid line displays the function, 1/ tanh(1.74y / (T c /T) — 1) 
w A /Abcs(^) = ■y F /A B cs(r)£o- In the dirty case (r n = 10 A , i.e., the 
mean free path I = vp/2T n = 0.05£o)> £i(T) behaves like the clean- limit 
BCS coherence length ~ v-p / Abcs(T') below T ~ 0.6T C , and is almost con- 
stant at low temperatures. The increase of the vortex core radius £i(T) with 
increasing T at high temperatures for the dirty case r n = 10 Ao, is more 
gradual than the temperature dependence of the clean-limit BCS coherence 
length (the solid line in Fig. 1(b)). This behavior is qualitatively consistent 




Fig. 1. The vortex core radius £i(T) (points) in the case of the s-wave 
vortex (|20|) as a function of the temperature T for several values of the 
impurity scattering rate r n . Lines are guides for the eye, except for the 
solid line in (b). In the plots, £1 and T are normalized by £o an d T c , 
respectively. Here, T c is the superconducting critical temperature, £o is 
defined as £o = ^f/Ao, and Ao is the BCS gap amplitude at zero tem- 
perature in the absence of impurities, (a) T n = 0, O.lAo, and Ao- (b) 
T n = Ao and 10 Ao- The solid line in (b) is a plot of the function, 
l/tanh(l. 7V (T c /T) - 1) « A /A B cs(T) = t> F /A B cs(T)£o, which repro- 
duces approximately the temperature dependence of the clean-limit BCS 
coherence length. 
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Fig. 2. The vortex core radius £i(T) (points) in the case of the chiral 
positive p-wave vortex ()2o|) as a function of the temperature T for several 
values of the impurity scattering rate T n . Lines are guides for the eye. In 
the plot, £i and T are normalized by £o and T c o, respectively. Here, T c q is 
the superconducting critical temperature in the absence of impurities, £o is 
defined as £o = vf/A-o, and Ao is the BCS gap amplitude at zero temperature 
in the absence of impurities. 



with a dirty limit result reported by Volodin et oZ.™ Note that the overall 
temperature dependence of £i in the dirty case is quantitatively different 
from that of Ref. 19 (£ e ff)> which is probably caused by the difference of 
the definitions of the vortex core radius. As displayed in Fig. 1(a) for the 
moderately clean case (r n = O.lAo, i.e., I = 5£o) and even in the relatively 
dirty case (r n = Ao, i.e., I = 0.5£o)> the vortex core radius £i(T) shrinks 
approximately linearly in T with moderate curvature below T ~ 0.6T C and 
saturates towards a finite value in zero-temperature limit. This gradual sat- 
uration due to impurities is in contrast to a sudden truncation of the KP 
effect which happens below a certain temperature related to the discrete 
energy levels of the low-lying vortex bound statesP^^ 



4.2. Chiral P-Wave Vortices 



In the chiral p-wave superconductors and generally in unconventional 
superconductors, the superconducting critical temperature decreases in the 
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r/rco 

Fig. 3. The vortex core radius £i(T) (points) in the case of the chiral 
negative p-wave vortex (|22|) as a function of the temperature T for several 
values of the impurity scattering rate T n . Lines are guides for the eye. In 
the plot, £1 and T are normalized by £o and T c o, respectively. Here, T c q is 
the superconducting critical temperature in the absence of impurities, £o is 
defined as £o = vf/A-o, and Ao is the BCS gap amplitude at zero temperature 
in the absence of impurities. 

presence of impurities. The units of the temperature T in Figs. 2 and 3 is 
T c q (the superconducting critical temperature in the absence of impurities). 
We obtain £i from the dominant components A:? (r) in Eqs. (|22[) and (|23|) . 

We show in Fig. 2 the vortex core radius £i(T) in the case of the positive 
vortex A( p ) [Eq. (|23|)]. At low temperatures, shrinks approximately 

linearly in T with moderate curvature and saturates towards a finite value 
in zero-temperature limit, as in the case of the s-wave vortex. The vortex 
core radius £i expands with the increase of r n , owing to the increase of the 
coherence length £ with the suppression of the pair potential far away from 
the vortex core. 

Figure 3 displays £i(T) for the negative vortex A( n ) [Eq. (|22j) ]. In con- 
trast to the cases of the s-wave vortex and the positive vortex, the vortex 
core radius £i(T) strongly decreases even in the presence of impurities, and 
shrinks toward almost zero for the moderately clean cases below r n ~ 0.2Ao- 

We can understand, in the following way, the reason why the vortex 
core shrinkage is robust against impurities. In the Eilenberger equations © 
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and the term which includes the impurity self energy is written as 

[a,g] = [E,g}=^-[(g),g], (44) 

where we referred to Eqs. © and (|12j) . Therefore, if the Fermi-surface- 
averaged Green function {g) is equivalent to the original Green function g 
(namely, if {g) = g), the above term (|44|) becomes zero and the impurity 
does not play a role in the Eilenberger equation. Such a special situation 
occurs in the negative vortex (|22|) of the chiral p-wave phase, and not in 
the positive vortex (|2HJ) or other usual vortices) 47 * 48 * It corresponds to a local 
restoration of the Anderson's theorem inside the vortex core IMIi $ l 4 5' l 4 7 l 4 $ l 4 v ) I W I 
This negative vortex (j22[) is more favorable energetically than the positive 
vortex ()23|) at least in the clean limitp^and therefore the negative vortex is 
likely to exist in the chiral p-wave superconductors. 



5. SUMMARY 

We investigated the temperature dependence of the vortex core radius 
£1 (T) defined in Eq. Q , incorporating the effect of nonmagnetic impurities 
in the Born limit. The isolated single vortex in the isotropic two-dimensional 
system was considered for the s-wave pairing symmetry and the chiral p-wave 
pairing symmetry. 

In the case of the s-wave vortex (|20|) . as seen in Fig. 1, at low temper- 
atures the vortex core radius £i increases with the increase of the impurity 
scattering rate T n up to r n ~ Ao owing to the saturation of the KP effect 
due to impurities. In contrast to it, at high temperatures the vortex core 
radius £i decreases with the increase of T n owing to the decrease of the co- 
herence length £ due to impurities. In the case of the s-wave vortex (|2()ll of 
the moderately clean state (r n = O.lAo, i.e., / = 5£o) and of the relatively 
dirty state (r n = Ao, i.e., I = 0.5£o)> as seen in Figs. 1(a), at low tempera- 
tures the vortex core radius £i(T) shrinks approximately linearly in T with 
moderate curvature and saturates towards a finite value in zero-temperature 
limit. This gradual saturation due to impurities is in contrast to a sudden 
truncation of the KP effect due to the discreteness in the energy spectrum of 
the low-lying vortex bound statesPoEDln the case of the chiral p-wave pair- 
ing system, as seen in Fig. 2, for the positive vortex (j23j) the shrinkage of the 
vortex core radius £i(T) saturates towards a finite value in zero-temperature 
limit by impurity scattering, analogous to the s-wave vortex. For the neg- 
ative vortex (|22lh however, the local restoration of the Anderson's theorem 
inside the vortex cor d 44 ' 45 ' 4 ^ ' 47 ' 48 ' 4 " ' 5 '^ yields a KP effect little affected by im- 
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purity scattering and, as seen in Fig. 3, the vortex core radius £1 (T) strongly 
shrinks linearly in T at low temperatures even in the presence of impurities. 

It would naturally be highly desirable to establish the KP effect exper- 
imentally beyond the present level. Our analysis shows that impurity scat- 
tering which is harmful for the KP effect in conventional superconductors, 
can under certain condition be harmless in a chiral p-wave superconductor, 
making this kind of system a good candidate for experimental tests. On the 
other hand, we expect that rather weak shrinkings of vortex core observed 
in NbSe 2 (Refs. MliMiWMty and CeRu 2 (Ref. E2J) may be explained partly 
in terms of the impurity effect, i.e., the vortex core shrinkage approximately 
linear in T with moderate curvature and its saturation towards a finite value 
in zero-temperature limit (Fig. 1(a)). Rather large extrapolated values of 
the vortex core radius at zero temperature observed experimentally might 
be partly attributed to effects of multiple gaps. Finally, we mention a multi- 
gap effect on the KP effect. Our preliminary results for ro(T) in a two-gap 
model show that a contribution from the Fermi surface with a smaller gap 
to the total supercurrent density j(r) around a vortex makes the position ro, 
at which \j(r)\ has its maximum value, shift outward away from the vortex 
center r = 0, leading to a finite ro at zero temperature in spite of the clean 
limit.^ Our detailed results for the KP effect in a two-gap superconduc- 
tor will be reported elsewhere. Sr2Ru04 (Refs. I43|ril|()2|) and NbSe2 (Ref. 
1^3*]) have multiple bands and may be effectively two-gap superconductors. 
MgE>2 is a typical two-gap superconductor.^ Further investigations on the 
vortex core shrinkage in terms of the multi-gap effects are left for future 
experimental and theoretical studies. 
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